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We study strong electron tunneling in the single-electron box, a small metallic island coupled to an 
electrode by a tunnel junction, by means of quantum Monte Carlo simulations. We obtain results, 
at arbitrary tunneling strength, for the free energy of this system and the average charge on the 
island as a function of an external bias voltage. In much of the parameter range an extrapolation to 
the ground state is possible. Our results for the effective charging energy for strong tunneling are 
compared to earlier - in part controversial - theoretical predictions and Monte Carlo simulations. 



CD 



I 

O 

o 



> 



O 

OO 

0^ 



I 

O 

o 



I. INTRODUCTION 

In the last few years charging effects in low capac- 
itance tunnel junctions have been studied extensively, 
on one hand because of their possible practical appli- 
cations in single-electron devices, on the other hand 
because of interesting theoretical questions related to 
them. These mesoscopic tunnel junctions are generic ex- 
amples of macroscopic quantum systems with discrete 
charge states and dissipation They can be fabri- 

cated by modern lithographic techniques with junction 
capacitances as low as ~ 10""'^^ - 10~^^ F. In these 
systems single-electron tunneling (SET) is strongly in- 
fluenced by the Coulomb interaction at temperatures 
T < Ec/kB ~ 1 - 10 K. The simplest device displaying 
these effects is the so-called single-electron box, where a 
metallic island is connected via a capacitor Cg and a tun- 
nel junction with capacitance Cj to a gate voltage source 
Vg @| • This external voltage polarizes the system, which 
allows a continuous tuning of the properties of the box. 
The (bare) charging energy of the system is given by a set 
of parabolas E{nG, n) — Ec{nG~nY as a function of the 
integer number, n, of excess electrons in the island and 
the continuous dimensionless gate charge hg = CGVc/e- 
The energy scale and curvature at small ug is given by 
Eg = e^/2C. It depends on the capacitance of the island 
C = Cg + Cj, and e, the electronic charge. The charging 
energy is minimized if n is as close to ug as possible. 
Hence at low temperatures, fcsT <^ Eg, and very weak 
tunneling in high resistance junctions, Rt S> Rk = h/e^ , 
n increases in a stepwise fashion as rtc is increased. In 
the range —1/2 < ug < 1/2, which for definiteness we 
consider in the following, the ground state corresponds 
to n = 0, i.e. Eo{nG) = E{nG,0)- At finite tempera- 
tures higher charge states are excited, and the expecta- 
tion value {njT has smeared steps, approaching a linear 
riG-dependence at high T. 

In the weak tunneling limit one can proceed pertur- 



batively in an expansion in the dimensionless junction 
conductance 
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In lowest order, transport by sequential single-electron 
tunneling processes is described by the 'orthodox the- 
ory' On the other hand, the perturbation 
theory also shows that equilibrium properties, such 
as the ground-state energy Eo{nG) and the low- 
temperature expectation value of the charge (n) = hg — 
l/{2Ec) dEQ{nG) / duG get renormalized by electron tun- 
neling processes. The latter is observable as a weakening 
of the Coulomb blockade ||J^ In first order in at one 
finds for the two quantities pl^pdl 
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and (for -1/2 < nc < 1/2) 
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We observe that even for weak at the corrections are most 
pronounced near the degeneracy points ug = ±1/2. In 
fact, at the degeneracy points and low temperatures any 
order perturbation theory fails. 

In order to cover stronger tunneling and/or a closer 
vicinity of the degeneracy points various theoretical ap- 
proaches have been pursued. Systematic higher order 
expansions in at have been performed for equilibrium 
|l2| , ^ , as well as transport properties ||lj]. These ex- 
pansions describe the system well for weak to inter- 
mediate tunneling strength as long as the degeneracy 
points are avoided or the temperature is not too low. 
Renormalization-group (RG) techniques [ p^JTs] as well 
as a partial summation of arbitrary order processes, ac- 
counting for resonant tunneling phenomena pM, cover 
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the regime of stronger tunneling. However, these ap- 
proaches have concentrated on two adjacent charge states 
and require specifying a high frequency cut-off of the or- 
der of the charging energy, thus introducing an uncer- 
tainty which prohibits a quantitative comparison. This 
hmitation has been overcome in the recent work of Konig 
and Schoeller within a 'real-time RG' approach. By 
including higher charge states they obtain cutoff inde- 
pendent results. Concentrating on the renormalization 
of the density matrix they can further cover nonequi- 
librium and transport properties. A cutoff-independent 
expression for EQ[nQ) had also been obtained within a 
'non-crossing' approximation scheme 

When Rt is lower, closer to Rk, the electron tunnel- 
ing leads to strong fluctuations of the charge on the is- 
land for all values of uq. In the limit of high tunneling 
conductances, 1/at can be treated as a small parameter 
and nonperturbative effects emerge. Also in this regime 
a RG analysis has been formulated Instanton tech- 
niques have been developed for the equilibrium proper- 
tics as well as a real-time saddle point expan- 
sion for transport properties ||2^. All these approaches 
p9|-p3[ arrive at the conclusion that in the strong tun- 
neling regime the charging energy E'^ is renormalized as 



E*c = f{at)Ec eM^2TT^ at) 



(4) 



While there exist general consensus about the exponen- 
tial dependence of E^ on at, much controversy remains 
about the pre-exponential factor. The instanton analy- 
sis of Ref. yields f{at) oc at ioi E^ < ksT <^ Ec 
(which is confirmed by the quasiclassical approach |20|) 
and fiat) oc a1 in the limit ksT <C E^. In contrast, in 
Ref. |22|] a cubic dependence f{at) oc af has been found, 
while Ref. reports a linear dependence for all tem- 
peratures. This controversy is one of the motivations for 
us to revisit the problem. 

The expression demonstrates that values of at ^ 
0.1 already correspond to strong tunneling, resulting in 
a substantial renormalization. (Accordingly, a dimen- 
sionless conductance with different numerical coefficients, 
more suitable for the strong tunneling regime could be 
defined. To avoid confusion, we prefer to use the weak- 
tunneling expansion parameter, as defined in Eq. (|^), 
throughout this paper.) In principle Coulomb blockade 
effects survive at T = for any value of at- For instance, 
clear signs of Coulomb blockade have been observed for 
at ^ 0.84 0. But the experimental observation of such 
effects for even larger at requires exponentially low tem- 
peratures ksT < E^, and therefore is hard to achieve. 
Thus, it is of most practical interest to investigate elec- 
tron tunneling for values of at of order 1 and less. For 
this purpose, and in order to fill the gap between the lim- 
its covered by the analytic approaches, we analyze the 
problem numerically by Monte Carlo (MC) techniques. 

MC results for the renormalized charging energy E^ 
had been obtained before by different groups p5| , ^ , p^ , 
and in the limit of small at there is a very good agree- 



ment among the recent data |23||2J] and those of per- 
turbation theory [Q. However, deviations arise for 
at ^ 0.2. They increase with increasing at and have 
a systematic character, hard to explain within the error 
bars. One possible source for the discrepancy for large 
at ^ 0.4 could be the fact that the simulations of Ref. 

were done down to substantially lower temperatures 
(typically keT > 2 x lO^^^:^) than those of Ref. |2|] 
(/csT > lO^'^-Ec)- But this does not account for the 
differences at smaller at ^ 0.4 where the temperature 
ksT ~ lO"^^^^ was found to be sufficient for conver- 
gence. The published MC data are also not sufficient 
to resolve the discrepancy concerning the prefactor in 
Eq. d). 

Other physical quantities of interest are the ground 
state energy of the system Eo{nG) (i.e. the lowest energy 
band) and the average charge on the island, (n), as a 
function of the gate charge for general values of the 
tunneling strength at- These quantities have been stud- 
ied analytically in the limits of weak [ p^|jTl| , [l^ and very 
strong tunneling jljjl^, at finite (not too low) temper- 
atures psi, 26|, as well as by the real-time RG analysis 
of Ref. |17| which appears to cover a larger parameter 
range. It is clearly of interest to extend these investiga- 
tions to intermediate values of at and low temperature, 
and to control the analytic results by numerical means. 

This paper is devoted to a detailed MC analysis of the 
single-electron box in the regimes of weak to strong tun- 
neling, < at < 1. In Sec. II, we describe the model and 
the computational method employed in our calculations. 
The MC results for the renormalized charging energy at 
low and finite temperatures are presented in Sec. III. In 
Sec. IV we present the free energy F{nQ) as a function 
of the gate voltage, and we show results for the mean 
charge on the island. A discussion of our results and of 
their relation to other work is given in Sec. V. 



II. MODEL AND COMPUTATIONAL METHOD 

A. Basic formalism 

The grand partition function of the single-electron box 
can be represented in terms of the path integral 

/oo 
dipo exp(27rimnG) x 

m— — oo 

/ 2?(^exp(-5M) . (5) 
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Here the "phase" variable ^{t) is conjugate to the island 
charge, and m is the "winding" number of the compact 
variable ip. The effective action S{i^\ is given by 
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where the first and the second terms account for the 
charging energy and the electron tunneling, respectively. 
The kernel a(T) reads in Fourier representation 



0-{t) = ~^ X! I^nl exp(iru;„) , 



(7) 



for Matsubara frequencies smaller than the electronic 
bandwidth, |wn| < D. For D > fcsT, which will 
be considered here, the kernel takes the form air) = 
at {irkBTy / ain^ {nkBTr). Apart from the bare charging 
energy scale Ec, the dimensionless conductance at, and 
the temperature, the model depends on the gate charge 
TT-G, which is proportional to the applied gate voltage. 

The effective action and partition function can be 
rewritten in terms of the phase fluctuations ^(t) = 
(P{t) — 27rmT//3, with boundary condition 9{0) = 0{(3), 
in the form 



Z = 



E 



(8) 



The coefficients Im{at,(3) = J VO exp {—Sra[0]) are to be 
evaluated with the effective action Sm[6{T)] = S[9{t) + 
2'KmT / [3]. They depend on the winding number m, the 
temperature, and the tunneling conductance, but are in- 
dependent of the gate charge uq- Thus, from a com- 
putational point of view, the problem reduces to calcu- 
lating the relative values of /m(at,/3), which can be ob- 
tained from the Monte Carlo simulations apart from an 
overall normalization constant. Since the partition func- 
tion is even and periodic with respect to nc. Zinc) = 
Z{—71q) = Z{nQ + 1), we can restrict our discussion to 
the region < < 0.5. 

For low tunneling conductance, the island charge is 
quantized and at zero temperature, the average number 
of excess electrons in the box, {n) , is a staircase function 
of the external voltage. In general, at finite temperatures 
and arbitrary tunneling the average charge can be deter- 
mined from the free energy F = —ksThiZ of the island 

by 
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From the gate- voltage dependence of the free energy one 
can, further, define a temperature-dependent effective 
charging energy for the single-electron box as 



E*c{T) 
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The zero-temperature limit of this quantity, E^ = 
Eq(0), is the renormalized charging energy discussed in 
the introduction. It coincides with the classical energy 



scale Ec for weak tunneling {at ^ 1), but is renor- 
malized in systems with stronger tunneling conductance. 
At high temperatures, the free energy F{nG) depends 
weakly on no, and the curvature E^{T) approaches zero. 
By using Eqs. (^) and (|To|), this effective charging energy 
Eq{T) can also be expressed as 

E*c{T) = 27r2fcBr(m2)„^=o , (H) 
where {m'^)na=o is a moment of the coefficients Im{at,l3). 



B. Monte Carlo method 

MC simulations have been carried out by the stan- 
dard discretization of the quantum paths into N (Trot- 
ter number) imaginary-time slices p7[ |. In order to keep 
roughly the same precision in the calculated quantities, 
as the temperature is reduced, the number of time-slices 
N has to increase as 1/T. We have found that a value 
N — A^Ec is sufficient to reach convergence of Jm, even 
for the lowest studied temperatures and for strong tun- 
neling, where convergence has been reported to be slower 
p^ . Thus, the imaginary-time step employed in the dis- 
cretization of the paths is Ar = /3/7V l/(4i?c)- This 
means that the high-energy cut-off oJc associated with 
this discretization is oJc ^ 27r/AT w 2^Ec- Repeating 
the calculation for a few data points with higher values 
of N and high-energy cut-off did not change our results. 

The classical Metropolis MC sampling |Q has been 
used to obtain finite-temperature results, and to extrapo- 
late, where possible, to T = 0. The partition function has 
been sampled according to Eq. (||) for temperatures down 
to ksT = i^c/SOO. These low temperatures were neces- 
sary to determine the zero-temperature effective charging 
energy for large tunneling conductance (see below). A 
simulation run proceeds via successive MC steps (MCS). 
In each MCS, all path-coordinates are updated. At each 
studied temperature, the maximum distance allowed for 
random moves was fixed in order to obtain an acceptance 
ratio of about 50%. For each set of parameters (at, T), 
we generated ~ 3 x 10^ quantum paths for the calculation 
of ensemble-averaged values. The starting configuration 
for the MC runs was taken after system equilibration 
at the considered temperature. In general, equilibration 
runs of about 2 x 10^ MCS were sufficient, but in some ex- 
treme cases, especially for strong tunneling conductance, 
equilibration runs of about 1 x 10^ MCS were necessary. 

In some cases, the acceptance ratio for jumps be- 
tween different winding numbers during a MC run be- 
comes very low, and sampling by direct jumps is ineffi- 
cient. This happens, in particular, at high temperatures 
{ksT > Ec), where Im/Io is very small for m 7^ 0, and at 
low temperatures (especially for strong tunneling), where 
the most relevant phase paths for different winding num- 
bers are very different. In these cases, we have calculated 
I-m' by carrying out simulations for fixed winding num- 
ber (say m), and evaluating in this fixed-m ensemble the 
average value 
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(12) 



for m =^ m' . The average value Rmm' so defined coincides 
with the ratio Im' / Im , since by definition we have 



Ijr^ ^ J VP cxp i-SM exp - S^, [9]) 



(13) 



By definition we have Rm'm = 1/ Rmm' ■ We have checked 
that the MC simulations satisfy this consistency relation 
within the statistical noise. We have also checked for 
some sets of parameters {T, at) that this method of cal- 
culating relative values of /,„ gives the same results as the 
direct Metropolis method, in which the winding number 
changes during a MC run. 

The MC method allows us to calculate the partition 
function Z as a function of T, at, and no- There is, 
however, a limitation on the parameter range that can 
be studied by this method. As noted in earlier publica- 
tions |l^,Q, due to the term exp{2Trimnc) in Eq. (||), 
the ratio Z{nc)/Z{0) for 7^ rapidly approaches 0, 
as the temperature is lowered. Since the MC method 
provides the coefficients with a limited accuracy, the 
partition function Zino) can be determined in a reliable 
way only when the ratio Z(nQ) / Z{Q) is larger than the 
numerical error in /„. The temperature range where re- 
liable results can be obtained for a given value of uq ^ 
reaches to lower values as at is increased, since the larger 
is at, the slower Z{nc)/Z{Q) decreases at low T. We 
found empirically that the lowest temperature for which 
the whole range < < 0.5 can be studied by the 
present method with sufficient accuracy, scales roughly 
as ksTmin ~ Eq{{))/20. On the other hand, this limi- 
tation does not apply for nc ~ 0. Hence, the effective 
charging energy E^{T) can be evaluated down to much 
lower temperatures. 



III. EFFECTIVE CHARGING ENERGY 

A. High-temperature regime 

At high temperatures {ksT 3> Ec), the partition func- 
tion of the single-electron box, including the effect of elec- 
tron tunneling, can be approximated by an expression 
similar to the classical result [E2| 



Z ~ ^ exp 

n— — 00 
00 



m— — 00 



-f3Ecin-nG)^ 
2-KimnQ 



PEc{T) 



(14) 



The temperature-dependent parameter Ec (T) takes into 
account quantum fluctuations. It reduces to the bare 
charging energy Ec both in the weak tunneling limit and 
at high temperatures. In the second form we made use 




0.1 0.2 0.3 0.4 0.5 

Inverse temperature (PEq) 

FIG. 1. The parameter Ec{T) is plotted in the high-T re- 
gion, for several values of the tunneling conductance at- Sym- 
bols indicate results of the Monte Carlo simulations: at =0.1 
(circles), 0.3 (squares), 0.5 (triangles), and 1 (diamonds). 
Error bars are smaller than the symbol size. Lines corre- 
spond to Eq. (|l^) , derived from semiclassical calculations in 
Refs. SEi. 



of a Poisson resummation to establish the relation to the 
winding number representation. At sufficiently high tem- 
peratures, where only winding numbers to = 0, ±1 have a 
non-neghgible contribution {Im/h ^ Ii/Io for > 1), 
Ec{T) is well defined. In this case Wang and Grabert 
p2t expressed the result of a semiclassical calculation, 
where Gaussian fluctuations around the classical paths 
are allowed, in the form (M). They find 



Ec{T) 



Ec 



1 2 at (3E, 



c 



(15) 



As the temperature is lowered, Ec{T) decreases as quan- 
tum fluctuations become more prominent. At still lower 
temperatures it is not guaranteed that the partition func- 
tion can be parameterized in the form dlj). 

If the partition function is of the form (llj) and Ec [T) 
can be defined, the latter is related to the effective charg- 
ing energy E*c{T), Eq. (|I§, by 



E, 



47r^ 

c{T) = — exp 



PEc{T) 



(16) 



This implies that at high temperatures, E^{T) vanishes 
proportional to Texp[— 7r^/(/3£'c')]- 

We determined Ec{T) by fitting the coefficients Im 
obtained from the MC simulation to a Gaussian profile. 
At sufficiently high temperatures, where Im converges 
quickly as a function of m, this is obviously possible. For 
this range the results are shown in Fig. 1. Our data points 
(symbols) are compared with the expression (|l5|) (contin- 
uous lines) for several values of at- An increase in the 
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0.01 0.1 

Temperature (kgT/EQ) 

FIG. 2. Effective charging energy Eq{T) as a function 
of temperature for several values of the tunneling conduc- 
tance. From top to bottom: at = (continuous line), 0.1 
(squares), 0.2 (open circles), 0.3 (black circles), and 0.4 (dia- 
monds). Dashed lines are guides to the eye. Error bars are of 
the order of the symbol size, unless shown explicitly. 



strength of tunneling and quantum fluctuations leads to a 
decrease of Ec- We observe good agreement between MC 
and analytic results up to (3Ec ~ 0.2, even for the high- 
est junction conductance considered. For larger values 
of /3 the high-temperature approximation ( p^ becomes 
insufficient, dropping below the MC results. 

In the temperature range shown in Fig. 1 the temper- 
ature dependence of the effective charging energy i?^ (T) 
is dominated by the exponential term in Eq. (p^), i.e. we 
find vanishingly low values for E'^^ irrespective of the 
strength of the tunneling conductance. 



B. Low-temperature regime 

We now turn to the low-temperature region, where the 
quantum fluctuations of the island charge dominate the 
thermal fluctuations and are responsible for the reduc- 
tion of the effective charging energy E^{T). The lat- 
ter can be conveniently calculated from the mean-square 
winding number obtained from MC simulations at finite 
temperatures, see Eq. (11). Our results are displayed 



in Fig. 2 for several values of the tunneling conductance 
at = 0, ...,0.4. In general, Eq{T) grows as the temper- 
ature is lowered, and reaches a plateau, the renormal- 
ized charging energy E^,{0), at low T. This saturation 
of Eq{T) should be found for any value of at- How- 
ever, for growing at lower and lower temperatures are 
required, consistent with the corresponding decrease of 
E^{0), which sets the scale for the convergence. 

At low temperature one can also employ an alternative 
method of calculation of Eq, and thus provide a consis- 
tency check. To do so, we first expand the free energy of 



the system in powers of uq, 



F{nG, T) = E*c(T) + Ai{T) n% + Ae{T) n% - 



(17) 



Taking a derivative of the partition function in Eq. (|^) 
with respect to we find 

MT) = ^ [3(™'>'-(^'))]„,.o ' (18) 

MT)^^[{m'')-15{m^){m')+30{mY]^^^^ , (19) 

and analogous expressions for the higher order terms. 
Provided that the low-tcmpcrature behavior is regular 
we can conclude from Eq. (jl^) that for (3E^{0) 3> 1, 
the average (m^) diverges as (m^) cx f3. Similarly, one 
has (m*) cx (3'^ and (m^) oc in this limit. Since the 
coefficients A„(r) are finite at all temperatures including 
T = 0, one obtains 

{m*)/{m^y = 3 + 0{l/f3), {m'^)/{mY = 15 + 0{l/(3), 

(20) 

and similar relations for the higher moments of m. This 
implies that at low T the distribution of the coefficients 
/,„ must be close to a Gaussian function of m, up to high 
winding numbers |m| < /3E^{0). In other words, in the 
limit of low T we can rewrite Im in the form 



Im oc exp[— a2(r) ~ 04 (T) — 



(21) 



with 02(r) - l/{l3Ec) and a4(T) - l/{f3Ec)^. Hence, 
we have ai{T) / a2{T) ~ l/{(3Ec) (the latter relation can 
be also derived from the fact that Im{T) scales at low 
T as J„(r) = I[m/^/PE^). Combining Eqs. (|) and 
(|l|) and defining Ec{T) = Tr^/l3a2{T) one obtains in the 
low-temperature limit 



(22) 

This relation is valid only for small values of nc, since for 
UG ~ 1/2 the sum converges slowly and higher values of 
TO, for which the Im are no longer Gaussians, gain impor- 
tance. Note that this expression for the partition function 
is formally identical to the high-temperature form given 
in Eq. (|l|). 

It is easy to show that in the limit T ^ the two 
quantities Ec{0) and E^{Q) coincide with each other. 
Indeed, in this limit the sum in Eq. ( p^ ) can be replaced 
by an integral, giving (m^)„g=o — f3Ec{0)/ {2n'^) and 
the identity Ec{Q) = Eq{0) becomes obvious. At low 
but finite T one finds from Eq. (02) 



EciT) 



l-ApEc{T)eM-PEc{T)] 



(23) 



5 



< 


' 


h 


0.5 

} 


□ 

I 




r 


t 


• 



0.001 0.01 0.1 

Temperature (kgT/E^,) 



FIG. 3. Convergence of the low-temperature effective 
charging energy Eq{T) as a function of the temperature, 
for a tunnehng conductance at = 0.5. The results of the 
present MC simulations are represented by black symbols. 
Open squares are results for Ec{T) derived from our Monte 
Carlo simulations. For comparison, results for Eq{T) found 
in earlier MC simulations for the same tunne ling conductance 
are given: Diamonds, data of Wang et al. open circle, 

data point of Hofstetter and Zwerger Error bars are 

shown when they are larger than the symbol size. The hor- 
izontal dashed line indicates the value to which our results 
converge at low temperature. 

The relations described above are well reproduced by 
our numerical MC analysis. Consistent with the relation 
(p^), at low T the coefficients Im are well described by 
a single Gaussian function, up to high values of m, char- 
acterized by only one parameter Ec{T). This parameter 
is presented in Fig. 3 as a function of temperature for 
at ~ 0.5. Also shown are the results for E^{T). The 
data points for Ec(T) are indicated by open squares, 
those for the effective charging energy E^{T) by black 
circles. Both converge to each other on a temperature 
scale set by Ec{T), consistent with (p3|). We also note 
that both functions converge to a common plateau at low 
temperatures, which defines E^{Q). Similar convergence 
is found for all other values of at; the results for E^{T) 
and Ec{T) derived from our MC simulations are indistin- 
guishable (within the statistical noise) for temperatures 
lower than ~ i;^(0)/(5fcB). 

Clearly, the two methods are not independent from 
each other. They are equivalent once the Gaussian dis- 
tribution of the coefficient Im has been established. How- 
ever, the fact that at low temperature the data display 
the Gaussian distribution with the required accuracy pro- 
vides a consistency check for our MC procedure. The 
combination of both approaches increases the reliability 
by reducing the chance of systematic errors. 

For comparison we also present in Fig. 3 the MC re- 
sults of two different groups for the same tunneling con- 




Tunneling conductance (a,) 

FIG. 4. Low-temperature effective charging energy, 
Eq{Q): (a) in the conductance region up to at — 0.4, and (b) 
for strong tunnehng. Black symbols: results of the present 
MC simulations. Data from earlier simulations are repre- 
sented by open symbols: triangles, from Ref. and squares 
from Ref. [^. The dashed-dotted and continuous lines cor- 
respond to second and third order perturbation theory in at, 
respectively. The dashed line represents the real-time RG of 
Ref. jl7j . The bold and dotted lines in (b) correspond to non- 
perturbative calculations in Refs. [fl9| and [E2], respectively. 
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ductance: diamonds ||2J,^ and open circles All of 
them are compatible with those obtained here except for 
those found in Ref. Q for PEc = 500, which lie clearly 
higher than our general trend. 

The results obtained for E^{0) from the low- 
temperature limit of E^(T) and Ec(T) are shown in 
Fig. 4 as black symbols. For comparison, we also present 
analytic and earlier MC results. In Fig. 4(a), we show the 
data for low to intermediate values of the tunneling con- 
ductance, < a* < 0.4, and in part (b) we present the 
effective charging energy E^{0) for the strong-tunneling 
regime on a semilogarithmic plot. 

Our MC results (black symbols) follow closely those 
given by Hofstetter and Zwerger in Ref. (open tri- 
angles), while for at > 0.3 the MC data of Wang et al. 
|P4| (open squares) are systematically higher than those 
found here. For completeness we add that very recently 
Goppert et al. carried out MC simulations in the 
charge representation, which implies an expansion in at, 
rather than the phase representation employed here and 
in previous MC simulations [ p3|j2^ ]. Since no error bars 
were given in that work, with symbol sizes too large to 
allow a distinction of different results, we do no compare 
with our results at this stage. 

In the figure we also present the analytic results of per- 
turbation expansions in the tunneling conductance up to 
second (dashed-dotted line) and third order (continuous 
line). In third-order, Goppert et al. found for the 
zero-temperature renormalized charging energy 



En 



= 1 - 4 at ^2 a? + ^3 a? 



0{at) 



(24) 



with A2 = 5.066 and ^3 = -1.457. We also show 
the Qft-dependence of the renormalized charging energy 
£^(0), as obtained from the real-time RG calculations 
||17| (dashed line). Our numerical values for i?c(0) agree 
with the results of the perturbation expansion and the 
RG approach only for at ^ 0.25. 

In Fig. 4(b) we display results for E^[Q) for stronger 
tunneling and compare with the predictions of sev- 
eral analytical calculations. In this regime the weak- 
tunneling expansions fail, while the real-time RG re- 
sults [0 (dashed line) show a qualitatively correct trend 
throughout. Already for at ^ 0.5 our MC results get 
close (within the error bars) to the result obtained by a 
nonperturbative instanton calculation [p9[ (bold line) 



E*c{Q) = 167r'*a?£;c exp(-27r2at -t- 7) 



(25) 



where 7 is Euler's constant. Since the numerics con- 
firms the result ( p5|) already for at ~ 0.5, and since the 
accuracy of the instanton analysis |l9) should increase 
with at, we expect that ( p5| ) remains accurate also in 
the regime at > 0.6 not covered by our MC simula- 
tions. For comparison, we also present in Fig. 4(b) the 
result for Eq{0) found in Ref. j2^ for strong tunneling 
at r = (dotted line), which predicts a pre-exponential 
factor ^ at . It lies more than one order of magnitude 
higher than the results of our MC simulations. 




0.2 0.3 

Gate charge (Hq 




0.2 0.3 
Gate charge (n^ 

FIG. 5. Free energy vs. the dimensionless gate voltage, 
no as obtained from Monte Carlo simulations, (a) Data for 
junction conductance at = 0.3 and different temperatures. 
From top to bottom, ksT/Ec = 0.01, 0.03, 0.05, 0.1, and 0.2. 
The width of the free-energy band, F(nG=0.5) — F{nG=Q), 
decreases as the temperature increases, (b) Same at fixed 
temperature ksT = 0.03i?c, for several values of the tunnel- 
ing conductance. From top to bottom: at = 0.5, 0.3, 0.1, 
and (dashed line). For each at, the free energy has been 
normalized by its value at no = 0.5. 



IV. GATE-CHARGE DEPENDENCE 

A. Free-energy band 

The MC method employed here allows us to deter- 
mine the "free-energy band" F{nG) as a function of the 
gate voltage, no, at finite temperatures. As an example, 
we display in Fig. 5(a) results at different temperatures 
for a dimensionless conductance at — 0.3. The error 
bars of the Monte Carlo simulations associated with the 
free energy are largest around uq = 0.5. At high tem- 
peratures, F{nQ) approaches a cosine shape, F{nQ) ~ 
E^lT)[l - cos(27rnG)]/(27r^), as expected for the clas- 
sical limit. At low temperatures, the free-energy band 
is closer to a parabolic shape, which is what we expect 
in limit of vanishing junction conductance at <C 1, with 
deviations which are most pronounced near = 1/2. 
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shown by the dashed line in Fig. 6. The comparison with 
our Monte Carlo results reveals differences. While the 



FIG. 

Width of the free-energy band, W = F(nG=0.5) — F{nG—0), 
as a function of the junction conductance at at several tem- 
peratures. Open symbols are finite-temperature results for 
kBT/Ec = 0.4 (squares), 0.2 (triangles), 0.1 (circles), and 
0.03 (diamonds). Black symbols correspond to the lowest 
temperature reached in our MC simulations. The continu- 
ous line is a fit to the expression W = Wo exp(— J^Qt), with 
K = 7.6 ± 0.2. The dashed line is the prediction for the 
low-temperature band- width for strong-tunneling ||l9|. 



A complementary picture is provided in Fig. 5(b). Here 
F{nc) is plotted for fixed temperature, ksT = 0.03Ec, 
while at is varied. The band is close to a cosine function 
for large conductance, and clearly differs from this shape 
for small at- For the ease of comparison in this figure the 
free energy has been normalized for each at to its value 
at no — 0.5. 

The temperature dependence of the bandwidth 
W{T) = F(nG=0.5) - F{nG=0) is shown in Fig. 6 for 
several values of at. For fcsT > Ec, one has W{T) ~ 
E^{T)/tt^ ~ 4kBTexp[-TT^/{l3Ec)], as follows from the 



high-temperature expression for E^,{T) given in Eq. (16) 
In general, one finds that W{T) increases as the temper- 
ature decreases, converging to a low-temperature limit. 
But this convergence is slower than that found for the 
renormalized charging energy E^{0). The black sym- 
bols in Fig. 6 correspond to the lowest temperature we 
could study for each at value. Within the tempera- 
ture range analyzed here we reached saturation of W{T) 
only for at ^ 0.3. In this range the low-temperature 
results can be fitted to an exponential law of the form 
WiT = 0) = Woexp{~Kat) where Wq = -Bc/4 and 
K = 7.6 ± 0.2 (solid line). This is still consistent, within 
the error bars, with the expression derived for weak tun- 
nehng W = Wo{l — Satin 2). For strong tunnehng, 
one expects for the bandwidth a similar scaling as a func- 
tion of at as for the renormalized charging energy, i.e. 
W{T = 0) ~ exp(— 27r^at). This form, with prefactors 
as predicted by the instanton calculation of Ref. , is 



exphcit result of Ref. |19| is not supported, the expo- 
nential dependence may be reached at larger values of at 
than covered so far. 

Note that even for weak junction conductance (where 
a good convergence of W{T) to its low temperature value 
is demonstrated), the second derivative F"{nG) close to 
no = 0.5 is still changing at the lowest temperatures con- 
sidered here. This is in line with the results of various 
analytic calculations [^- ^^ |l6 , 1^] , where a logarithmic 
divergence for F"(nG=0.5) was found at T = 0. While 
this divergence cannot be seen directly in our finite- 
temperature MC simulations, its precursor is clearly ob- 
served: F"{nG—0-5) increases continuously as tempera- 
ture decreases for all at values considered here. 



B. The charge on the island 

The average number of excess electrons in the island, 
(n), is of practical interest, since it can be measured di- 
rectly by measuring the voltage of the box. It follows 
from the free energy by Eq. M). In Fig. 7(a) we display 
(n) vs. the gate charge no at several temperatures for 
a dimensionless conductance at = 0.1. As expected, at 
high temperature the average charge follows closely the 
line (n) — uq, while it decreases for lower T as Coulomb 
blockade effects become more pronounced. According to 
the numerical restriction discussed above for the calcu- 
lation of the partition function Zinc), the lowest tem- 
perature that could be studied for all values of no at 
at = 0.1 was ksT ~ O.OSEc- For these parameters the 
average charge {niriG)) has reached the zero-temperature 
limit for < uq ^ 0.35, as can be seen from the con- 
vergence of the curves in the figure. For gate charges 
no closer to 0.5, no saturation is found for the temper- 
atures covered, and the T = value will be lower than 
our finite-temperature results. Our conclusions are sup- 
ported by the real-time RG calculations by Konig and 
Schoeller @- They find a zero-temperature {n^nc)) 
close to our low-temperature results up to no ~ 0.35. 
However, consistent with earlier work |10-l^Jl^, the RG 
calculations predict at T = a logarithmic divergence 
for the slope of {n{na)) at no — 0.5. This is beyond the 
range covered by the MC data. We note, however, that 
the lowest temperature presented in Fig. 7(a) is not far 
from the lowest temperatures presently attainable in the 
laboratory, as for typical values of Ec/ks = 1 K and T 

20 mK, one has fcsT = 0.02Ec. 

In Fig. 7(b) we present the dependence (n) vs. uq for 
several values of the tunneling conductance, as obtained 
from MC simulations at a temperature T = 0.03Ec/kB- 
In this figure, the tunneling conductance increases from 
bottom to top, and for the largest at shown {at = 0.4) 
the dependence is close to a linear one. For compari- 
son, we give also the average charge for the case at = 
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0.1 0.2 0.3 0.4 0.5 

Gate charge (n^) 




Gate ctiarge (Hq) 

FIG. 7. Average number of excess electrons, (n), in the 
single-electron box. (a) For a dimensionless conductance 
at — 0.1 at various temperatures. From top to bottom: 
kBT/Ec = 0.6, 0.2, 0.1, 0.04, and 0.03. (b) At tempera- 
ture T = 0.03i5c/fcs for several values of ctt. From top to 
bottom: at — 0.4, 0.3, 0.2, and 0.1. For comparison, the 
average charge for = is also given (dashed line). 



(dashed line). These finite-temperature results compare 
well with those found recently from Monte Carlo sim- 
ulations in the charge representation by Goppert et al. 
|l3| at a lower temperature {(3Ec = 10^). The main dif- 
ference between our finite-temperature results and those 
reported in ref. fl^ is found again in the region close to 
nc = 0.5, where the above-mentioned logarithmic diver- 
gence for d{n)/dnG is expected to appear. 

In order to analyze how long Coulomb blockade are ob- 
servable it is convenient to study the derivative d{n) /dna 
as a function of uq. For at ^ this derivative has 
its maximum at no = 0.5, and it approaches a plateau 
around no = 0, the height of which is related to Eq{0) 
by Eq. (|l^) . This decrease is a measure for the strength 
and observability of charging effects. The results of our 
MC simulations are presented in Fig. 8 for tunneling con- 
ductances in the range 0.1 < at < 0.5, at tempera- 
ture fc^r = O.OSi^c (a* decreases from top to bottom). 
For convenience we normalized d{n)/dnQ by its value at 
ug — 0.5 for each at- We see that even for the highest 
conductance presented, at = 0.5, charging effects remain 
observable in the range of temperatures comparable to 
the corresponding E^{0)/kB- the decrease of d{n)/dnG 
at riG = as compared to the value for uq = 0.5 is ~ 3%. 
A similar conclusion can be reached for other values of 
at, i.e. charging effects should be observable (at least) at 
temperatures of the order of E^{0)/kB even in the limit 
of large junction conductances. These results are com- 
patible with recent measurements by Chouvaev et al. , 
where clear signs of the persistence of charging effects 
were observed in a sample with an effective conductance 
as large as at ~ 0.84, even at temperatures larger than 
the renormalized charging energy. 



V. DISCUSSION 




0.2 0.3 
Gate ctiarge (Hq) 

FIG. 8. Derivative of the mean charge in the island with 
respect to the gate voltage, d{n) / dna, for different values of 
at at fcflT = O.Q'iEc- From top to bottom: at — 0.5, 0.4, 
0.3, 0.2, 0.1. The dashed line corresponds to the absence of 
electron tunneling (at = 0). For each at, the corresponding 
curve has been normalized by the value at no = 0.5. 



The detailed MC analysis carried out in the present 
paper allows us to determine the low-temperature val- 
ues of the renormalized charging energy E^ for a single- 
electron box in the range of weak to strong tunneling 
< Oft ^ 0.6. Our MC data interpolate well between the 
perturbative result s 021 and those of a nonperturbative 
instanton analysis |19|] in the limits of low and high at, 
respectively. Our data demonstrate that the third order 
perturbation theory |l^ yield quantitatively correct val- 
ues for E^ for at ^ 0.25. On the other hand, we find 
quantitative deviations at higher values of at- 

A similar conclusion has to be drawn concerning the 
validity of the RG approach of Ref. . This approach is 
not equivalent to a direct perturbative expansion ^2 13 
in at since it allows for a partial summation of diagrams 
in all orders in at- On one hand, for large at it works 
better than the perturbative expansion capturing 
the qualitatively correct trend of the renormalized charg- 
ing energy Eq to decrease with increasing at up to high 
values. On the other hand, the neglect of higher-order 
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vertex corrections in the RG may lead to quantitative 
errors at large at- While the explicit validity range of 
the RG analysis is difficult to establish analytically, the 
corresponding information can be extracted from a com- 
parison to the MC data. This cornparison reveals differ- 
ences between the RG approach and the MC data 
for at > 0.3. 

For stronger tunneling we expect a crossover to an ex- 
ponential dependence of on at- This is clearly shown 
by our MC data already at at ^ 0.5, where it approaches 
the result of Ref . ||T^ . This crossover is not displayed by 
the perturbative expansions in powers of at , and from a 
general point of view it appears that it cannot be cap- 
tured within any finite-order expansion. For this reason 
also MC algorithms formulated in the charge represen- 
tation may be limited to not too large values of at- On 
the other hand, the real-time RG result has been fit- 
ted to the exponential dependence (|^) for 0.5 ^ at < 1, 
but it requires a substantially different pre-exponential 
function /(a*) as compared to that found within the in- 
stanton approach |l9[ . 

Additional support for both the analytical results ob- 
tained in Ref. | [T9[ | and our numerical data in the strong 
tunneling regime comes from an independent MC calcula- 
tion of Ref. . These authors analyzed numerically the 
correlator {cos{lp{t) — ip{0))) as a function of temperature 
at different values of at , and observed a sharp crossover 
in the behavior at a temperature T*{at). Since this tem- 
perature should be close to the renormalized charging 
energy T*{at) ~ E^{0)/kB, it is interesting to compare 
both quantities for different values of at- This compari- 
son, carried out in Ref. revealed a very good agree- 
ment of the data for T*{at) with the result of the 
instanton calculation ( |25| ) (see Fig. 1 in Ref. ||2^), which 
in turn agrees with our MC data in the strong tunneling 
regime. 

Our numerical data for E^, are in agreement with the 
previous MC data of Ref. |23| , extending them to larger 
at and to substantially lower values of the temperature. 
Although the accuracy of our MC analysis is higher as 
compared to that achieved in Ref. (the error bars are 
smaller and the number of points higher), it is satisfac- 
tory to observe that our data show the same trend as that 
found in Comparing our results to the MC data of 
Ref. we find good agreement for low at ^ 0.3, while 
for larger values of at their data are systematically higher, 
clearly showing a different trend with increasing at as 
compared to the one indicated by our data. There are 
several reasons why we believe that this discrepancy can- 
not be ascribed to insufficiently low temperatures used in 
our simulation: 

(i) For all studied values of at ^ 0.6, the calculated 
E^{T) increases smoothly as the temperature is lowered 
and it clearly saturates at T < E^{0)/5 (see Fig. 2). The 
temperature T = 2 x lQ~^Ec /ks is sufficient to observe 
saturation for all at up to ~ 0.6. Moreover, a systematic 
difference between our data and those of Ref. exists 
already at relatively low tunneling strengths 0.2 < at < 



0.4, where the saturation of Eq{T) was observed by all 
groups already at T > 10~^Ec, i.e. well above the lowest 
temperature employed in our simulations. 

(ii) The quantity Ec{T) converges to the same value 
E^,{0), showing saturation at approximately the same 
T < £'p(0)/5 as E^(T). This convergence is clearly seen 
in Fig. 3 for at — 0.5. The same behavior of Ec{T) 
was observed for all other values of at. Since the val- 
ues Eq{T) and Ec{T) were calculated independently, 
the chance of systematic errors is reduced. Furthermore, 
these two quantities monotonously converge to the same 
value from below and from above, respectively, thus pro- 
viding a lower and an upper bound for E^{0) at each 
T. According to (^3|) these bounds should merge at suffi- 
ciently low temperatures. This is exactly what our data 
demonstrate. 

(iii) The lowest temperature reached in our MC sim- 
ulations and in those of Refs. | p^ , p9t is the same. We 
observe (see e.g. Fig. 3 for at = 0.5) that at higher tem- 
peratures UbT/Ec ^ 0.01 the data of Wang |Q are still 
fully consistent both wi th P3| and with our data. How- 
ever, the data point [|^j29j|at a lower ksT/Ec ^ 0.002 
is by more than a factor of 2 higher than all other points. 
Neither was this result confirmed by our MC analysis, 
nor do we see any physical reason for such a rapid jump 
of E'^{T) at low T, where this quantity should already 
approach its zero-temperature value. 

Very recently one more numerical study of the renor- 
malized charging energy in the strong tunneling regime 
was performed in Ref. jlj], using a different MC algo- 
rithm. The centers of the symbols for the data points 
|l3| lie in-between our data and those of Ref. . Unfor- 
tunately, the actual error bars are not presented in Ref. 
|l3| and the symbols are so large that they cover both 
our data and those of Ref. |2^ , thus making a detailed 
comparison impossible at this stage. 

In summary, the Monte Carlo method employed here 
has been shown to be well suited to study charging ef- 
fects in the presence of an external voltage at not too 
low temperatures. An important advantage of the MC 
analysis is that it covers both perturbative and nonper- 
turbative regimes and allows us to describe a crossover 
between them. We conclude that the combination of the 
expansion (M), the instanton result (25) and the numer- 



ical data covers all values of at, thus providing complete 
information about the renormalized charging energy E^ 
at low temperatures. Our results also resolve the con- 
troversy between the strong tunneling theories [l9| ] and 
p2| in favor of the former. The temperature range for 
which results for the average charge in the island (n) can 
be obtained for < na < 0.5, covers most of the tem- 
perature region actually studied in the experiments. The 
Monte Carlo simulations indicate that charging effects 
are observable in the single-electron box even for strong 
tunnehng, at temperatures of the order of E^{0)/kB- 
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